Theory of neutral and charged exciton scattering with electrons in semiconductor 

quantum wells 
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Electron scattering on both neutral (X) and charged excitons in quantum wells is studied 

theoretically. A microscopic model is presented, taking into account both elastic and dissociating 
scattering. The model is based on calculating the exciton-electron direct and exchange interaction 
matrix elements, from which we derive the exciton scattering rates. We find that for an electron 
density of lO^cm"^ in a GaAs QW at T = 5K, the X~ linewidth due to electron scattering is roughly 
twice as large as that of the neutral exciton. This reflects both the X~ larger interaction matrix 
elements compared with those of X, and their different dependence on the transferred momentum. 
Calculated reflection spectra can then be obtained by considering the three electronic excitations 
of the system, namely, the heavy-hole and light-hole IS neutral excitons, and the heavy-hole IS 
charged exciton, with the appropriate oscillator strengths. 

PACS numbers: 34.80.-i, 78.67.Do, 73.21. Fg 
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INTRODUCTION 



The broadening of exciton emission lines due to 
electran-exciton collisions was first observed in bulk 
GaAsld. An experimental comparison between exciton- 
electron ai 
bulk GaAi 



exciton-exciton scattering mechanisms 
and in GaAs quantum wells (QW's 
showed that the exciton-electron scattering efficiency is 
an order of magnitude larger than that of the exciton- 
exciton process. Both scattering processes are en- 
hanced for the 2D excitons as compared to bulk ex- 
citons. Theoretical calculations of elastic and inelas- 
tic exciton-electron scattering, were reported for bulk 
semiconductorscl and for QW'sB. The semiclassical treat- 
ment of the latter neglected the exchange term which was 
shown recently to be the dominant term in QW'sQ. 

The existence of the negatively charged exciton A"" 
(trion) was first proposed by LamperilS as the semicon- 
ductor analogue of the hydrogen ion. Due to their con- 
finement in the growth direction, 2D trions have a bind- 
ing energy which is an. order-of-magnitude larger than 
that of the bulk trionsLL-This fact facilit|ated the obser- 
vation of trions in CdTal3 and in GaAstll QW's contain- 
ing a low-density two-dimensional electron gas (2DEG). 
Theoretical calculations of the binding energy of both 
negatively {X~) and positively {X~^) charged excitons 
were performed using various trial pfLmrtions, showing 
that only the singlet state is boundO'EjuJ. Particular at- 
tention was devoted to the modification of the trion's 
properties in the presence of strong magnetic fields. In 
this limit, the X~ binding energy is increased, and the 
electron spin triplet state becomes bound as welliil. ft 
was also observed that the dependence on the magnetic 
field of the X~ and Ar+ binding energies, as well as their 
Zeeman splitting, differ drastically (while being nearly 
identical at zero magnetic field)lla. 

While these aspects of the trions where extensively 
studied, there are very few reports on their broadening 



mechanisms00. In this paper we present a theoreti- 
cal model for the scattering of both neutral and charged 
excitons with electrons in QW's. This paper contin- 
ues our previous work which introduced neutral excitons 
elastic scattering with electrons in the context of cav- 
ity polaritona3. Here we extend our model, to incorpo- 
rate both excitons and trions, taking into account elastic 
and inelastic scattering processes. The model is based 
on calculating the exciton-electron direct and exchange 
interaction matrix elements, from which we derive the 
exciton scattering rates. These are integrated over all 
final excitonic states, resulting in the exciton linewidth 
due to scattering, as a function of its initial momentum 
(or energy). A major difference between the exciton and 
trion scattering lies in the charge of the latter, which re- 
sults in a divergence of its matrix elements in the limit 
of zero transferred momentum. This divergence, origi- 
nating from the infinite range of the Coulomb potential, 
is treated by applying the Lindhard model for the po- 
tential screening. The screening action complicates the 
trion linewidth dependence on the electron density; thus, 
while for electron densities larger than 5 • lO^cm"^ the 
trion linewidth due to electron scattering in a GaAs QW 
at r = 5K is comparable to that of the neutral exci- 
ton, for a very dilute 2DEG (rie « 5 • lO^cm"^) the trion 
linewidth becomes roughly an order-of-magnitude larger 
than that of the neutral exciton. Inelastic scattering is 
also much more efficient in the trion's case, due to its 
smaller binding energy. The calculated reflection spectra 
are obtained by considering the three electronic excita- 
tions of the system, namely, the heavy-hole and light-hole 
IS neutral excitons, and the heavy-hole IS charged ex- 
citon (trion), with the appropriate oscillator strengths. 
A qualitative validation of our calculations is given by 
considering photoluminesence (PL) measurements that 
were done on a mixed type I - type II GaAs/AlAs QW 
(MTQW) structurellZl. 

The paper is organized as follows. In section II, we 
present in detail the model for the neutral exciton - elec- 
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tron scattering. In section II] we consider the trion- 



electron scattering. First we construct the trion's wave 
function, which is used in the matrix elements calcu- 
lations. We then treat the divergences that originate 
from the infinite range of the Coulomb potential, by ap- 
plying the Lindhard model for the potential screening. 
The effect of the screening is discussed in the context 
of the trion linewidth dependence on the electron den- 
sity. In section IV we use our model to calculate QW 
reflection spectra and discuss their relevance to the avail- 
able experimental data. A summary and conclusions are 
given in section In Appendices A and B, we pro- 
vide some details of the direct and exchange integrals 
calculations for the neutral exciton - electron scattering. 
Finally, a derivation of the trion binding energy, using a 
Chandrasekhar-type trial function, is given in Appendix 
C. 



II. A MICROSCOPIC MODEL FOR EXCITON 
ELECTRON SCATTERING 



fermionic Hilbert space of electron-hple pairs, using the 
notations of Tassone and YamamotoE£l. It is a superpo- 
sition of wave functions with different electron momenta 
ki and electron and hole z coordinates, given by 



ki 

(III) 



where cj^ ^ (dj^ is the electron (hole) creation op- 
erator with in-plane momentum kx and z^. (z^) coordi- 
nate, mo, Mx are the electron and exciton in-plane effec- 
tive masses, respectively, a = ma/M^, and (ji'k.i.^e, Zh) = 

A/" [l -I- (Afc)^] '^^'^ Xe{ze)Xh{zh) is the in-plane Fourier 
transform of the exciton wave function (see Appendix A) . 
A state comprising an exciton and an unbound electron 
having kc and will be written as 



In this section we present a detailed description of the 
exciton-electron scattering, considering separately the 
elastic process, for which the exciton remains bound, and 
the inelastic process, where the exciton breaks into an 
unbound electron-hole pair. 



|kx;ko) = ^ J dZedZhdZc4>*a]^^+i^^{Ze, Zh)lpl^{Zc) X 
ki 

cLk„..dL+k„.,4.,.jO), (11.2) 



A. Elastic scattering 

We first construct the state of a single exciton with 
an in-plane center-of-mass (CM) momentum kx in the 



where '(/'ko (-Zc) is the electron wave function. Applying 
the Coulomb interaction operators Vee, Veh to this state, 
we have for the electron-hole and electron-electron inter- 
actions: 



|kx; kc 



jkx! kc) 



= J2 J '^^'=^^'»'^^c^qi(^'» ~ {Ck,+ki+qi(2^e,Z,,)V'k_^(Zc)cLk,,^^4,+ki,^;,c[„,2, 
ki.qi 

+ Ck.+ki(^e,2/i)^ke(^c)cLki,..4.+ki-qi,z.4.+qi...} |0)- (II-3a) 
dzedzhdzcVq.ize-^ Zc)(t>li,^+i,^{ze,Zh)tp^^{zc)cl^^^_^^^,ji^+^^ (II.3b) 

1. „ 



ki.qi 



The first term in Eq. (II. 3a) represents the Coulomb interaction between the constituents of the exciton, thu s con - 
tributing to its self energy, and can be discarded in the calculation of the scattering matrix elements. In Eqs. ( II. 3a ) 
and ( II.3b| ), V^i(z) = ±^|^^e~*i^ is the 2-dimensional Fourier transform of the Coulomb interaction, where A is 
the QW area, and the plus (minus) sign is used in Vee (Veh)- Using the anticommutation relations for the fermion 
operators, one can easily find the scattering matrix elements to be 

(kx +q;kc - q|V;/l|kx;kc) ^ ^ C^^e^^/lC^Zc [V'q(zh - zj(?!)„(k,+q)+k^ (Ze, Zft)(/)*k^+ki (^e. ^^)V'ke--q(2c)'0k„ (^c)- 

ki 

Vk,+-k^-q{zh - 2c)0ak,+ki-/3q(2c, ^/OCk^+q-k^ (^e> Zh)llJk^-q(Ze)lPuS^c)] (II.4a) 
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kx +q;kc - q|V;e|kx;kc) = ^ / rf^erf^^/i^^^c [Vq(Ze - Zc)(/)ak,-/3q+ki (^^e, (^e, ^ft)V'ke-q(2:c)'0ke(^c)- 

ki 



Vki+ke-q(Ze - Zc)0ak,-,3q+ki(^:c,Z/i)Ck,+ki(2:e,2/OV'k,-q(2:e)V'k„(^c)] , (II.4b) 



where we denote [3 = I — a. The first term in each of 
Eqs. (II. 4) contributes to the direct (classical) Coulomb 
interaction and the second contributes to the exchange 
matrix element. The direct term is calculated analyti- 
cally (the details can be found in Appendix A) , resulting 
in: 

^dir('?) - HM) [5(A/?9/2) - g{Xaq/2)] , (II.5) 

where L is the QW width, A is a variational parameter 
associated with the exciton Bohr radius in the QW, and 
we have defined the dimensionless functions 



9iM) 
h{Xq) 



21-3/2 



[1 + iM) 

-1 + Lq 



(Ag)3 [1 + (Lg/2^)2 



(11.6) 



It is evident from Eq. ( |II.5| ) that the direct term is iden- 
tically zero for equal electron and hole masses (a = /3 = 
0.5). Similarly, we have for the exchange term 



yoxc(Ak,q) = 
[g(A|ki+q- 



Se^A^ 



d2fci.g(A|ki 



Ak|)-5(A|q-Ak|)]/i(Afci) 



aq — Ak|) x 
(11.7) 



where we have defined Ak = ko — akx. The exchange 



term (II. 7) is computed numerically, however its angular 
part can be calculated analytically (see Appendix B). It 
is convenient to transform to dimensionless direct and 
exchange integrals given by 



^ ^ 2 e^A3 



The direct and exchange integrals are plotted in Fig. |i| 
as a function of the transferred momentum g, for the 
case Ak = 0, where the angular dependence of /oxc dis- 
appears. As g — > the direct integral approaches zero 
, while the exchange integral has its maxinu«n (this is 
also the case for exciton-exciton interactionta) . In the 
general case, the exchange integral is a function of the 
transferred momentum q, the momentum difference Ak 
(which can be regarded, for convenience, as the in-plane 
momentum of the colliding electron in the rest frame of 
the exciton), and of the angle 6 = Z(q, Ak). The ex- 
change interaction term has the following features: (i) 
The interaction favors the case q = Ak. Physically, this 
means that the electron is inclined to transfer as much 
momentum as possible to the exciton, preferably in the 
same direction, (ii) The interaction retains its strength 
for quite large values of q (or Afc) even though the ex- 
citonic wave function vanishes much more rapidly with 



momentum, (iii) The differential cross section is largest 
for 9 = and decreases to a minimum for back scattering 
(9 = tt). 
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FIG. 1: Calculated direct and exchange integrals vs. trans- 
ferred momentum. The inset shows — /dir on an expanded 
scale. 



Next we calculate the scattering rate of an exciton, 
with an initial in-plane momentum kx, to an exci- 
tonic state with an in-plane momentum kx -I- q. The 
exchange matrix element (II. 7) depends on five vari- 
ables: kx, kc,q, cos 6, cos 4>, where we have defined (fi = 
Z(kx,kx -l- q). It is convenient to transform from the 
angles fl— and (j) to the angles 7 = Z(kx,q) and S = 



cosO = 



cos ( 



fco cos S — a/cx cos 7 



^kg + a^k^ — 2akxke cos(7 ± 6) 
kx + q cos 7 



2kxq cos 7 



, (Afc^O) 
, (kx^-q). (n.8) 



It can be verified that for the special case of fcx = 0, we 
have 9 = 5 and </> = 7. Conservation of energy and mo- 
mentum for the exciton-electron scattering process reads 

Ex{\kx + q|) + £;c(|ke - q|) = Ex{kx) + E,{k,), (II.9) 

where (Ec) is the exciton (electron) kinetic energy. 
This equation is satisfied with: 



1 



2Mx cos S 



[q{Mx + Too) + 2fcxmo cos 7] , (II. 10) 



The scattering rate is then found using Fermi's golden 
rule 
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The scattering rate is then found using Fermi's golden rule 



Wx-olkx kx + q) 



A 
nh 

An 



^2^ l^dir + K=xc|m-fce 



d£c(|kc-q|) _ dgo(fce 
dk(2 dkc 



/fd(fce) [1 - /fd(|ke - 



ka = k^ 



7rft3 |c0S(5| q 



^* iKiir + Kxci' /fd(fc:) [1 - /fd(ik: 



(II.U) 



where /fd(fcc) = [e(^<=(''''=)-^)/'^^3T ^ i jg ^-j^g electron 
Fermi-Dirac distribution function. In the two dimen- 
sional case that we are considering, the chemical poten- 
tial is ^ — kBT\og{e^f/^'^'^ — 1), and we assume that 
the Fermi energy is given by the free 2D electrons value: 
Ef = Trti^nc/me. We note that in the high temperature 
limit, Ef <^ kgT, the Fermi-Dirac distribution is prac- 
tically classical, and that for low electron densities, the 
chemical potential beco mes n egative 
the 



(II., 



the integrand in ( 11.11 ) 
function of the variables: 



Using the relations 
can be expressed as a 
/cx, 9, cos 7, cos 5. The exciton linewidth due to electron 
scattering is calculated by integrating over all final exci- 
tonic states 



= (fcx) 



hA 



qdqd-yw^^cik^, q, cos j) . (11.12) 



Figures || a,b show the heavy-hole exciton (M^ — 
0.177me) and light-hole exciton (Mx — 0.306me) 
linewidths at T = 807^ and T = 5K for = 5- lO^cm-^. 
The large linewidths obtained for relatively low ric reflect 
the high efficiency of the electron scattering mechanism. 
This should be compared to an exciton linewidth of w 
0.2meV for acoustic phonon scattering at T = 8QK. The 
larger linewidth obtained at T = 5K (Fig. ^d) is ex- 
plained by noticing that the exciton-electron interaction 
matrix elements favor small energy transfer transitions. 
Thus, at high temperatures, the electrons have too high 
energy to be effective scatterers. Decreasing T thus in- 
creases the scattering rate until the electron gas becomes 
degenerate, and fewer final states are available for the 
scattered electrons. In figure ||c the heavy hole exci- 
ton linewidth at fcx = is plotted as a function of ric 
for T = 5K and 80K (the light hole exciton linewidth 
exhibits a similar behaviour). It is seen that while at 
T — 80K the linewidth amplitudes scale practically lin- 
early with Tie for a large range of electron densities (up to 
no ~ lO^^cm^'^), the linearity region at T = 5K is much 
smaller. This linearity threshold seems to be in accord 
with the onset of the phase-space filling effect which be- 
comes noticeable at higher densities and effectively en- 
larges the exciton Bohr radiustj. As the temperature 
decreases, the effect of phase-space filling becomes im- 
portant at a much lower electron density. We note that 
the functional dependence of the linewidths on the ini- 
tial exciton momentum is almost unchanged within the 
linear regime. Increasing ric further results in a shift of 
the maximum linewidth from /cx = to higher momenta. 
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FIG. 2: (a) Heavy exciton (Mx = 0.177me) and light exciton 
(Mx = 0.306me) linewidths due to electron scattering for a 
2DEG with rio = 5 • lO^cm"^, as a function of the exciton 
initial in-plane momentum, at T = 80K ; (b) Same as (a) at 
T = 5K ; (c) Linewidth of heavy hole exciton with initial mo- 
mentum fcx = 0, as a function of ric, for the two temperatures 
considered. 



B. Exciton Dissociation scattering 

We now consider the case when the electron-exciton 
scattering results in dissociation. This scattering process 
is depicted schematically in Fig. ^. The initial state of 
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FIG. 3: Electron-exciton inelastic scattering scheme 



the system is given again by Eq. (CL3) while the final 
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state is that of three free particles, symbohcally written 
as |ke — q; k'; kx + q — k'), where q is the momentum 
transferred from the free electron to the electron-hole 
pair, and k' is the second electron momentmn. Calcu- 
lating the Coulomb interaction matrix elements between 
the initial and final states of the system, we again obtain 
direct and exchange interaction terms: 

Kit (kx,q,k') = ^ (^)'^\(Ag) X 

[g{X\aK + q - k'|) - g(A|ak, - k'|)] (11.13a) 

Kt^^(kx,q,k',ke) = (^] MA|k'-ke|) X 



[g{X\aK + q - k'|) - giX\aK + q k,|)](II.13b) 



where g{Xq), h{Xq) are given in (|L6|). We note that con- 
trary to the elastic scattering case, here the direct and 
exchange matrix elements are comparable. In order to 
calculate the scattering rate we first introduce energy- 
momentum conservation 

E^{k^) + E,{k,) = ^b + Sh(|kx + q-k'|) + 

^e(A:') + ^c(|ke-q|), (11.14) 

where Eh is the exciton binding energy. Denoting the 
angle v — Z(kx, k') together with 7, S which were defined 
previously, this equation is satisfied with 

^^-2/3gcosn ^ ' +q +k + 

2a [k^q cos'y — k'k^ cosv — k'q cos{-j — z/)]} (11.15) 

The scattering rate is then calculated to be: 



I dis 

X — c 



(kx 



-q) = 



k'dk'dvdS kl 



\vi:±v^ /M(fc:)[i-/fd(ik; 
[i-/td(fc')], 



|cos(5| q 

-q|)]x 



11.16) 



where the -I- (-) sign between the direct and exchange 
matrix elements corresponds to the singlet (triplet) elec- 
tron spin configuration. For an unpolarized electron gas, 
an averaging over the two spin configurations must be 
performed, in order to take into account the contribu- 
tions of all the electrons to the exciton linewidth due to 
dissociating scattering: 



(fcx) 



1, 



idis,+ 



(fcx). 



(11.17) 



rx!!e^(^x) (rx!!e~(^x)) denotes the exciton linewidth con- 
tribution from the singlet (triplet) spin configuration, 
and both are calculated using Eq. ([1.12). We note that 
the spin configuration averaging is irrelevant in the elas- 
tic scattering case, since the direct term there is much 



smaller than the exchange term. This results in simi- 
lar contributions from both singlet and triplet configura- 
tions. 
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FIG. 4: Heavy exciton (Mx = 0.177me) linewidths due to 
dissociating scattering with electrons, as a function of exciton 
initial in-plane momentum, with nc = 5 ■ lO^cm"^, at (a) 
T = SO/C ; (b) T = bK. The insets show the maximum 
linewidth dependence on ric. 

Figure ^ shows the heavy hole exciton linewidth due 
to dissociating scattering at rio = 5 • lO^cm^^, and 
T = 80K. Ahhough the magnitude of r^i!^(fcx) is of the 
same order as the elastic scattering linewidth (compare 
to Fig. Ija), its functional dependence on the exciton in- 
plane momentum is very different. In particular, we note 
that the maximal linewidth is obtained at a very large 
momentum (fcx ~ 5A~^). This is due to the fact that in 
order for an exciton with initially small fcx to be ionized, 
it must scatter on an electron with energy large enough 
to overcome its binding energy. This is less likely as the 
temperature decreases, as is evident from Fig. where 
T^l^ik^) is plotted at T = 5K. The insets in Figure | 
show the maximum values of rxifc(fcx) vs. Ue- As in the 
elastic scattering case, the linear dependence on Ue holds 
for much larger electron densities at T = 80K. 



III. CHARGED EXCITON (TRION)-ELECTRON 
SCATTERING 



At low temperatures and low ric, charged excitons are 
formed. In this section we calculate the trion linewidth 
due to scattering with free electrons. Three scattering 
processes are considered: elastic scattering, where the 
trion remains bound, capturing scattering, where the 
heavy-hole exciton captures a free electron to form a 
trion, and dissociating scattering where the trion dissoci- 
ates into a heavy-hole exciton and an extra free electron. 
The latter process is important due to the small bind- 
ing energy of the trion with respect to the exciton. We 
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neglect the scattering in which both electrons become un- 
bound. In order to facilitate the matrix elements calcula- 
tion-jye consider a simple two parameter Chandrasekhar 
typecj trial function for the trion 

Xh{zh) [0(rih)0'(r2h) ± 0'(rih)0(r2h)] , (ULl) 

where the + (-) sign a pplies to the singlet (triplet) spin 
configuration. In Eq. ( III.l ), nh — |ri — rh|, {i = 1,2) 
are the in-plane coordinates of the two electrons with re- 
spect to the hole in-plane coordinate, Zi are the electrons 
z coordinates, and Xe{zi), Xh{zh) are the confinement 
functions in the z direc tion, taken to be the same for all 
constituents (Eq. ([A-3|) ). The two electron orbitals are 
given by the normalized functions 



0'(r) 



ttA^ 



rA'2 



-r/A 



-r/A' 



(III.2a) 
(III.2b) 



results are fairly accurate, giving us the confidence to 
proceed in calculating the trion-electron scattering ma- 
trix elements using the wave function (III.l). 

We construct the state of a single trion with an in- 
plane C M m omentum kj, similarly to that of a single 
exciton ( II.l ): 



ki.ka 

Sl,S2 



Cs{su^2)cX^^,c%^J^^^^,+^^,jQ), (III.5) 
where we have added the spin index to the electron 
creation operators, and denoted at — mc/Mtr, and 
is{si, S2) = (•S'Isi, S2) as the projection of a generic spin 
configuration of two electrons on the singlet spin con- 
figuration. The in-pla ne Fourier transformed trion wave 
function in Eq. ( 111.5 ) is given by 



•^ki ,k2 (^1 , 22 , ^h) = MtrXe {zi )Xe {z2)Xh (^/i) 0ki 



ki (Pk2 ' 

(III.6) 



where, 



and the trion wavefunction normalization factor is given 

by 



Mr 



3/2 



(III.3) 



A 



1 



A [l + (Afc)2]3/2 

87rA'2 1 



[1 + (A'A:)- 



13/2 • 



(III.7a) 
(III.7b) 



where 



4AA' 



(A + A')2' 



(111.4) 



The trial function (III.l) was used inE3 for the case of 
a two dimensional negative-donor center (i.e., nih ~* 00 
limit). In addition, these authors performed the bind- 
ing energy calculations using a more complicated ver- 
sion of (III.l) which includes a correlation term (with an 
additional variational parameter). The variational pa- 
rameters A, A' are calculated by maximizing the trion 
binding energy for the two spin configurations. The 
details of the binding energy calculation are given in 
Appendix C, resulting in El'' = 0.985meV,A = 143A, 
A' = 300A(£;f = 0.765fieY,X = 15lA, A' = 7650A) 
for the singlet (triplet) spin configuration, where E"^"^ is 
calculated with respect to the heavy hole IS neutral ex- 
citon energy. We note that the triplet configuration is 
barely bound, in accordance with the experimental ob- 
servations, thus it can be disregarded. Comparing. with 
the much more elaborate treatment of Stebe et a£j, our 



A state comprising a trion and a free electron having ke, 
Ze, and Se will be written as 



|kt; ke) = 



E 

ki.ka 

Si,S2,Se 



dzidz2dziidzeQ{si, S2) x 



itr* 

'^atkt+ki,Q:tkt+k2 



(zi, Z2, Zh)lpt{Ze) X {lll.i 



ki, 21 C_"k2, 22*^14 +ki+k2, 2ft Ckjz JO)- 



In the following we shall use these wave functions to cal- 
culate the scattering matrix elements for the various scat- 
tering processes. 



A. Elastic scattering 

Applying the Coulomb interaction operato rs Ve i, Ve2, 
and Veh to the initial state of the system ( III.8| ), and 
discarding the self energy contributions, we can write 
a generic expression for the three scattering matrix el- 
ements: 
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kt + q; kc - q\Ve^\kt■, ko) = ^ ^1 ^ J dZtdz[dZhdz',^Vq,iZe - 2»)'/>a,(kt+q)+ki,at(k,+q)+k;, (^0 X 



kl,k2 '<1>'<2,S 



k'^ ,k2,qi ^1 '"^2' 

JJJ\ 



CsK,4)V'k,-q(4)Ck.+k,,..k.+k.(^^)G(si,S2)V'k.(^e)(0|ciC2C3cl44|0), - l,2,h), (111.9) 



where we put the 1/2 prefactor to indicate averaging over 
the initial electron spin states, and used a numeric index 
for the electron operators, indicating both spatial and 
spin degrees of freedom. Note that we have omitted the 
hole operators since they anticommute with the electron 
operators and do not contribute additional constraints. 
Using the electrons anticommutation relations, the Fermi 
vacuum expectation value of the operators reads: 



|ciC2C3444|0) 



'534('525'5l6 - (^15^26) + 



<524(^15<536 - Ssb^ie) + 314(635626 - 625636), (111.10) 



where we have symbolically written 6i 



3k,, k 



6(zi 



Zj)6g. g^. For the 1st electron interaction term we have: 

l = (ke-q;z^;s^) 2=(-k^;z^;s^) 

3 = (-k[;z[;s[) 4 = (-ki - qi; zi; si) (III.U) 

5 = (-k2; 22; S2) 6 = (-ke + qi; Ze; Se) 



The second electron interaction term is found by making 
the substitutions: ki — > ki — qi, k2 ^ k2 + qi in Eqs. 
( III. 111 ). Similarly, the electron- hole interaction term is 
found by substituting: ki ki— qi in (III. 11). Applying 
the delta functions of ( [II.IOD to Eq. ( [II. 9[ ), it is evident 
that the contributions of the z-dependent parts are iden- 
tical for all three terms, since we have taken identical 
con finem ent functions for the electron and the hole (see 
Eq. A- 3 ) . Assuming that the trion remains in the singlet 
electron spin configuration, we sum over the spin degrees 
of freedom, resulting in selection rules for the three elec- 
tron final spin states, which determine the relevant signs 
of the various terms. Performing the in-plane momentum 
integrations, and collecting the matrix elements into one 
direct and two exchange contributions, we finally find: 



Vdir(q) = Ksh(\q) { [g(Xatq/2)g(\'Ptq/2) + g(X' atq/2)g(\Ptq/2) + 2K,^g(\atq)g(\f3tq) 
g(Xatq/2)g(X'atq/2) + K^g^Xatq)] } 



and 

viic^(q,Aki 



(III.12) 



J d^kih(Xki) { [xX'Kg(Xatq)g(X\ki - Akt + atq\) + X'^g(Xatq/2)g(X'\ki - Akt + atq\) 



[5(A'|ki - Akt + q|) - g(X'\Akt ~ q|)] + [A ^ A'] x [A ^ A']} 



(111.13a) 



Vg}(q,Akt) = J d^hh(^h){[xX'Kg(X\k,-atq\)g(X'\q-Akt\) + X''g(X'\ki-atq\/2)g(X\q^Akt\) 



g(X\ki - Akt + atq\) + [X ^ A'] x [A ^ A']} 



(III.13b) 



where we have denoted f3t = l — ct t, Akt = k^ — atkt, and 
A = XX'/(X + X')^. In Eqs. (|III.12| )-( pLl^ ), the prefactor, 
multiplying the dimensionless quantities, is given by 



1 



1 ine^X^ 
■f k2 L^Aeo ■ 



(III.14) 



We note that the direct term ([11.12) diverges at q 0. 
This divergence originates from the infinite range of the 
Coulomb potential, which is manifest whenever a scat- 
tering event between two charged particles occurs (note 
that for the neutral exciton this divergence cancelled out. 



due to the equal contributions with opposite signs from 
the electron and hole constituting the exciton). In prac- 
tice the Coulomb potential is screened by the presence 
of the 2DEG. Since we are considering a low density 
2DEG, its Fermi wavevector, fcf — ^/2TTn^ is small, thus 
the semiclassical (Thomas-Fermi) approximation, which 
assumes that q fcf, is inadequatecfl. In particular, the 
two-dimensional semiclassical analysis results (at T ~ 0) 
in a density- independent screening wavevector qo , imply- 
ing that the presence of a single electron is sufficient to 
screen the external potential, which is clearly unrealistic. 
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We adopt instead the Lindhard approach which utihzes a 
perturbative scheme to evaluate the induced charge den- 
sity in first order of the total potential (By total we mean 
the sum of the external potential, i.e., that of the trion, 
and the 2DEG potential that is induced by the trion's 
presence). The in-plane Fourier transformed dielectric 
functionjjjjhich results from the effect of the screening is 
given bynil 



£s{q) = 1 + ^gs{q)as{q) 



where 



90 



2mee^ _ 2 



(III.15) 



(III.16) 



and flB is the electro ns effec tive bulk Bohr radius. Two 
factors appear in Eq. (III.IE), which influence the screen- 
ing. The first is the screening form factor given by 



-q\z-z'\ 



(III.17) 



which arises from the finite QW width. In the infinite 
barrier limit we are considering , th e confinement func- 
tions, x{z), are given by Eq. (A_-3), and gs{q) can be 
calculated analytically, resulting in 



9s{q) 



i-^i -1 + Lq 



5 



{Lqf 



3 

32? 



(Lq)' 



(qLy 



(III.18) 



This screening form factor is smaller than 1, thus it 
reduces the effect of the screening, compared with the 
strictly two-dimensional case {gsj L ) = 1). The sec- 



ond factor which appears in Eq. (111.15) is given by 
1 



asiq) 



2, /fd(k- 



/fd(fc) 



(k + q)2 - fc2 



(III. 19) 



We note that as{q) does not depend on the particular de- 
tails of the QW, but only on the two-dimensional nature 
of the 2DEG motion. Substituting the screening dielec- 
tric function into the Coulomb interaction, we find the 
screened potential to be 



1 27re2 



27re2 



es{q) Aeo q 



Aeo q + qo9siq)as{q) 

(III.20) 



thus the use of the screened potential implies replacing 
h{\q) in Eq. (|III.12|) with 



hs{Xq) 



h{Xq) 



1 



L'-'a-B 



hi\q)as{q) 



(III.21) 



In general, the effect of the quasi two-dimensional screen- 
ing saturates in larger electron densities, in contradistinc- 
tion to the bulk case, resulting in a complicated depen- 
dence of the trion linewidth on nr.. We note that since 



there are no divergences in the ex change terms, replac- 
ing h(Xki) with /is(Afci) in Eqs. (III. 13) amounts to a 
very small effect (less than 1% in the relevant parameter 
range) . 
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12 3 4 

FIG. 5: Calculated dimensionless matrix elements at T = 
5K for the case of Akt = 0, vs. transferred momentum: (a) 
Direct (rie = lO^cm"'^) and exchange terms (note that since 
the exchange terms are not affected by the screening, their 
dependence on rio is negligible) ; (b) Total matrix element for 
various electron densities ric . The inset provides a close-up at 
small q. 



Figure shows the direct and exchange dimensionless 
matrix elements as a function of transferred momentum 
q, for the case Akt = 0, where the angular dependence of 
Vcxc , Vcxc disappears. The effect of the screening on the 
total dimensionless matrix element is shown in Fig. ^jb, 
where various electron densities are considered, giving 
rise to a change in the screening action. As expected, 
the effect of the screening is noticeable in the region of 
small momenta of the direct matrix element. 




10 12 



FIG. 6: Calculated trion linewidth due to elastic scattering vs. 
initial trion momentum, at T = 5K for two electron densities. 



9 



The elastic scattering rate of the trion, wt-o(kt 
kt + q) , is calculated by writing the energy-momentum 
conservation equa ti on an d finding its solution fc*, simi- 
larly to Eqs. (n.£)-( II.f 1 ). Integrating over all final tri- 
onic states, we calculate the trion linewidth due to elas- 
tic electron scattering, rt-c(fct), plotted for two electron 
densities in Fig. ^. The linewidth dependence on ric is 
quite complicated. First, as we are considering T = 
the density range for which the linewidth exhibits a linear 
dependence on rig is quite small, similar to the neutral 
exciton scattering (cf. Fig. ||c). Second, the screening po- 
tential becomes larger as the density decreases, thus for 
very small densities (rie ^ 5 • 10^) decreasing ric results 
in a larger trion linewidth. We note that for such low 
densities, it is probable that other screening mechanisms 
become appreciable (e.g. lattice impurities), thus such a 
behavior will probably not be observed experimentally. 



B. Electron capture 

We now consider the case in which a heavy-hole exciton 
captures a free electron thereby forming a trion. Th e ini - 
tial and final states of the system are given by Eqs. (11.2) 
and ([II. 5), respectively. The calculation of the Coulomb 



interaction matrix elements between these states essen- 
tially follows the procedure described in the previous sec- 
tion. In this process there is no exchange term. Instead, 
the two electrons which constitute the trion form a sin- 
glet spin configuration. Defining Ak = atkt — axkx, the 
calculation yields: 



Vcap(Afc) ^AAcap / d^fci/l(Aofci) 



XX' 



(Ao + A')2 



g (A'|ki + Ak|) -g(A'Afc) 
AA' 



5 A 



ki 



Ak 



(Ao + A)2 
<?fA' 



g^Ajki + Ak|) -g(AA/c) X (III.22) 



Ak 

ki + 



where we have defined the reduced Bohr radii: A = 
AoA/(Ao -I- A) and A' = AqA'/ (Aq -I- A'), and the prefactor: 



^A 



(III.23) 



Fig. ^ shows the dependence of the dimensionless ma- 
trix element on Afc. Proceeding as before, we write con- 
servation of energy: 



£;x(fcx) + ^e(kt-kx) = ^t(kt 
which is satisfied for 
1 



(III.24) 



k: 



fcx cos 6 ± \ 13, 



Elm, 



klsrn^e] (III.25) 



where 9 = Z(kx,kt). 
serving solutions, fc*, 



We note that for the energy con- 
the interaction matrix element is 



constant since Ak(fcj*) = ^ axCttEl^nic/'h? . Thus, the de- 
pendence of the heavy hole exciton linewidth on its initial 
momentum, shown in Fig. |^b, comes solely from energy 
momentum conservation and the Fermi-Dirac distribu- 
tion of the scattered electron: /fd(ke) = /fd(kt — kx). 
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FIG. 7: (a) Calculated dimensionless matrix element of the 
capturing process vs. Afc at T = SA"; (b) X^h calculated 
linewidth due to capturing scattering for rtc — lO^cm"^ at 
three temperatures. 



C. Charged exciton dissociation 

Lastly, we consider the case in which the scattered 
trion dissociates into a neutral exciton and a free electron. 
This scattering process is shown schematically in Fig. |^. 
The initial state of the system is given by Eq. ( [II. 8| ) while 
the final state is that of three free particles, symbolically 
written as |kx; kc; kt-l-q — kx), where q is the momentum 
transferred from the free electron to the trion, and kx is 
the neutral exciton CM momentum. 



(kt,Cs(si,S2)) 


(kx, Sx) 






t •^^^ 






(kt + q-kx,s') 




e • «- 


e • 






e • ^^^^ 


(ke , ) 


(ke - q, Se) 



FIG. 8; Electron-trion inelastic scattering scheme 

Proceeding as before, we calculate the interaction ma- 
trix elements between the initial and final states of the 
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system. To this end we denote the final spin state of the 
three electrons as |sx,s',Sg) (see Fig. ||). Assuming the 
interaction is spin independent, this final state will be one 
of the two states: (s, m) — =^5)1 where the state with 
m=^(m = — i) corresponds to the initial electron spin 
Se being up (down). Actually, we need not concern our- 
selves with the exact form of the final electron spin state, 
since the projection of the initial spin state (which is a 
tensor product of a singlet state with the additional free 



electron spin) on a generic final spin state will extract the 
relevant contributions. The proper (anti)symmetrization 
of this final spin state is performed inherently by the an- 
tico mmutat ion relations of the electron operators, as in 
Eq. ( p^II.lC ). As in the elastic scattering case, we find 
one direct and two exchange contributions. Denoting 
Akt = ko — atkt, and Akx = ko — axkx as before, these 
terms are expressed as: 



Vdir(Ak,q) 
and 

y«(q,Akt,Akx) 



K£?(q,Akt,Akx) = 



XX' 



XX' 



(Ao + A)2 



g{X'^^gCXAk)+g{X' 



Ak 



ttx 



.g(A|Ak + q|)-g(AAA:) 



(Ao + x'y 



A ^ A' 
A^ A' 



(III.26) 



AA 



-AAdi^y d^hh(Xoki)g{Xo\ki - Akxl) 



5 A 



ki- 



Ak 




Ak 






q 


— + q 


)] 


X 


Qfx 


Qfx 







<?(A'|q - Aktl) + g(A|ki - Akt + q|)<? A' 





Ak 




— + q 


)} 




Qfx 





--^^fdishs { Ao 



Akx Akt 



at 



XX' 



(Ao + xy 



+ {A^ A'} 



ff(A'|Akt-q|)3 A 



(III.27a) 



Akt - /3xAkx 



g(AAfc) 



5(A'— )-5(A'|Akt 



AA' r A ^ A' 
(Ao + A')2 1 A ^ A' 



q 

(III.27b) 



r 



where A, A' and Ak were defined in the previous sub- 
section and the pr efactor multip lying the dimensionless 
quantities in Eqs. ([II.26)-([II.27), is given by 



16V2e^Xt 

A/dis — 



2tt 



3/2 



(III.28) 



The d irect term ([11.26) and the 2nd exchange term 
( |lll.27f| ) diverge at 9 ^ and at (Akx/ax - Akt /at + 
q) — > 0, respectively, implying the need of using a 
screened Coulomb potential, as was done i n the t rion 
elastic scattering. Note that the direct term ( III. 26 ) de- 
pends only on fct,fcx,g, and on the angles 7 = Z(kt,q) 
and v = Z(kt,kx), whereas the exchange terms depend 
also on the initial electron momentum fco and on the an- 
gle S = Z(kc,q). In practice, we express everything in 
terms of q, Afct, A/cx, and the angles: 61 — Z(Akt, Akx), 
6*2 = Z(Akx,q), and 6*3 = Z(Akt,q), given by: 



cos til ~ 

cos 02 — 
cos 6*3 = 



1 



(fc^ - a^k^kc cos(7 - 6 - v)- 



AfcxAfct 

atktkf, cos(7 — 6) + a-^atk^kt cos ly) (III. 29a) 
kc cos S — axkx cos{ j — v) 



Afcx 

fee cos (5 — Qftfct COS 7 



(III.29b) 
(III.29c) 



The calculation of the scattering rate is performed as 
before, where conservation of energy now reads: 

Et{kt) + E,{k,) = + Sx(fcx) + £;c(kt + q-kx)-H 



Sn(k, 



(III.30) 



Equation (111.30) is satisfied for 
1 



1 



gcosJ [ h? 
q [kt cos 7 - /cx cos(7 - i^)]} 



9' + ^ (/3tkt - kx)' 



(III.31) 




FIG. 9: Calculated trion linewidth due to inelastic scattering 
with electrons vs. its initial momentum, for ric = lO^cm"^. 



Using Fermi's golden rule as in Eq. (11.17), we find the 
scattering rate u't_o(kt ^ kt -I- q). Integrating over all 
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final trionic states, we find the trion finewidth due to 
inelastic scattering as a function of its initial momentum 
kt, r^'i\(kt), shown in Fig. | for = lO^cm-^. 



IV. EXCITONS AND TRIONS IN QW'S 
CONTAINING A 2DEG 



In this section we apply the linewidth calculations for 
the various scattering processes, presented in the previ- 
ous sections, in order to calculate the line shapes of neu- 
tral excitons and trions. These line shapes can be directly 
compared with measured refiection spectra, taken from 
MTQW structures containing a variable density 2DEG. 
Without electrons in the QW and for T < lOOK, the 
two dominant scattering processes are due to acoustic 
phonons and static disorder. We note that the radiative 
lifetime of a free QW exciton is Trad ~ 25ps, resulting 
in a radiative homogeneous broadening of 7rad ~ 26/ieV. 
The exciton - acoustic phonon interaction leads to a ho- 
mogeneous broadening whose Lorentzian line shape, Cq, 
is characterized by 7phon = 0.2meV for the Xhh and 
7phon = 0.28meV for the Xih a^-fT = 80K, using the 
parameters of GaAs/AlAs QW'sE3. The IS exciton en- 
ergy, Eis, has a Gaussian distribution (characterized by 
7dis) due to the roughness of the QW surface. Convolv- 
ing these-two distributions results in the familiar Voigt 
functioned 




1 2 
E(meV) 

FIG. 10: A schematic picture of the convolution resulting 
from electron scattering. The solid lines represent various 
exciton initial line shapes To{E — E'). The contributions 
of states with fcx > to the fcx = state are weighted by 
the value of the Lorentzian £x-c {E' , T^-c{E')), given by the 
dashed line. 



The calculated line shapes are identified with the imag- 
inary part of the dielectric function in the QW that 
is related to the exciton (trion) resonance. Use of the 
Kramers-Kronig relations yields the real part of the di- 
electric function. Examples of calculated exciton line 
shapes together with the associated real part of the di- 
electric function are depicted in Fig |ll| for two electron 
densities. 



ME) 



7ph, 



(27r)3/27di, {E - E^Y + (7phon/2)^ 



(IV.l) 



where we take i?is = 0. When a low density 2DEG is 
present, Tq{E) is convolved with a Lorentzian lineshape 
that is associated with the homogeneous broadening that 
is due to the electron scattering: 



I{E) 



ME 



E')C^^, {E',T^ME'))dE' , 

(IV.2) 



where rx_c(-E'') is given by Eq. ( 11.12 ), with the exciton 
energy, £", corresponding to k^ . A pictorial description 
of the convolution in Eq. ( [V.2[ ) is given in Fig. hq, for 
both elastic and dissociating exciton-electron scattering 
processes (the trion-electron scattering exhibits a similar 
behavior). It is seen that the electron scattering admixes 
excitons having fcx > with the k^ — Q state. The degree 
of admixture is determined by the value of the Lorentzian 
that is due to electron scattering, given in the figure by 
the dashed line. Note that since rx-c(-E'') is a decreasing 
function of E' , the Lorentzian peak is shifted to higher 
exciton energies. (This is observed experimentally as an 
increased shift with increasing nj3). In the calculation, 
the Gaussian width, associated with the inhomogeneity, 
was taken as a fitting parameter, yielding the value Udis = 
0.2meV for both Ahh and Xy^ (the same value was taken 
also for the trion line). 
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FIG. 11: (a) Heavy-hole (Mx = 0.177me) exciton normalized 
line shapes due to scattering with electrons for two electron 
densities; (b) The corresponding dielectric function real part. 

Combining the three exciton line shapes (heavy-hole 
and light- hole neutral excitons, and heavy- hole trion), by 
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using the appropriate oscillator strengths, the reflection 
spectra of a QW can be obtained. For high temperature 
(T = 80K), where the trion is dissociated, these calcu- 
lated spectra and their dependence on ric were shawn to 
be in good agreement with the experimental dataCI. 

In the low temperature regime {2K < T < 20K), 
the vindication of our results with experimental measure- 
ments is harder to perform. Since the heavy hole exciton 
and trion lines are very close in energy, and the trion 
oscillator strength is much lower than that of the exci- 
ton, it is practically impossible to discriminate between 
them by reflection measurements. PL experiments were 
performed in MTQW with 2DEG, at low temperatures, 
facilitating the observation of the trion lineEZl, but there 
the detailed balance between the exciton and the trion 
line intensities is determined by the 2DEG density, and 
the system dynamics should be addressed by solving, e.g., 
rate equations. 
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FIG. 12: Trion's linewidth vs. electron density, for three 
temperatures due to (a) elastic scattering ; (b) both elastic 
and dissociating scattering. 



In Fig. 



12fi the trion linewidth due to elastic electron 
scattering is shown as a function of no for three temper- 
atures. The combined effect of both elastic and dissoci- 
ating electron scattering on the trion linewidth in showu. 
in Fig. |r^. Qualitatively, the trion line observed intZl 
indeed broadens considerably at much lower than the 
exciton line, as indicated by our results. The electron 
densities stated inliU are probably an overestimate, as 
implied by the fact that the neutral exciton linewidth 
remains unchanged up to rio = 3 • lO'^^cm"^, in clear 
contradiction to both our calculations and other experi- 
mental data (e.g.|]). It should be noted that calibrating 
the photoexcitation intensity /l to in MTQWs is not 
an easy task. Even if one assumes a linear dependence of 
ne on (which is valid for a limited range), the slope de- 
pends on the temperature, since the holes tunneling times 
decrease exponentially with T, so a different calibration 
should be performed not only for different experiments, 
but also for each temperature. Finally, the strong en- 



hancement of the dissociating scattering with increasing 
temperature, evident in Fig. |l2|b, implies that the trion 
line should practically vanish for T > 10 K at very low 
electron densities. 



SUMMARY AND CONCLUSIONS 



In this paper we presented a theoretical study of the 
various exciton-electron scattering processes that take 
place in a QW which contains a low density 2DEG. We 
have demonstrated that even for small densities, of the 
order of rie ~ 5 • lO^cm"^ in the case of a GaAs/AlAs 
QW, exciton-electron scattering is much more efficient 
than both exciton-phonon and exciton-exciton scatter- 
ing, thus playing a central role in the system dynam- 
ics. The dependence of the exciton linewidth on ric was 
shown to be drastically different for the cases of a clas- 
sical (T = S,OK) and a degenerate (T = 5K) 2DEG. 
In general, the trion-electron scattering poses a much 
harder problem, one of the major difhculties being the 
evaluation of the screening effects. We showed that the 
screening model has a crucial effect on the trion-electron 
scattering matrix elements, and hence on the resulting 
trion linewidth. This is particularly true for very low 
densities {ric < 5 • lO^cm"^), where the 2DEG screen- 
ing is less effective, making other possible mechanisms 
more dominant. Using the calculated scattering rates we 
showed that dissociating scattering can produce a siz- 
able broadening of the excitonic lines, in particular for 
the trion, whose binding energy is considerably smaller 
than that of the neutral exciton. 

A simple method to produce excitonic line shapes from 
the calculated linewidths was devised, allowing our re- 
sults to be easily compared with experimental reflection 
spectra. A reliable source of experimental data is achiev- 
able from reflection measurements taken from a MTQW 
structure, embedded in a microcavity (MC). In the limit 
of large MC mode energy, the exciton linewidths are re- 
stored, since their coupling with the MC mode is barely 
noticeable. A major advantage of MC experiwients is that 
the trion line is easily resolved by reflectior£3, allowing a 
direct comparison with theory. 

The calculated interaction matrix elements can be 
readily used to evaluate polariton-electron scattering 
rates in MCs with 2DEG. In a recent experiment 
we have demonstrated polariton final-state stimulation, 
assisted by these scattering mechanisms, which was 
hitherto oba/pyed using polariton-polariton scattering 
mechanismEilcj. Incorporation of the various scattering 
processes which were considered in this paper, in a con- 
sistent explanation for these nonlinear effects, will be the 
topic of a forthcoming publication. 
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Appendix A 

In this appendix we provide the details of the direct 
term calculation. Separating the coordinates in the QW 
plane {x — y) from the perpendicular coordinate (z), and 
denoting the electron and hole in-plane momenta by ki 
and k2 , respectively, we write the in-plane Fourier trans- 
form of the exciton wave function $[(re|| i -^e), (r;i|| , Zh)] 

'^iiiMize.Zh) = ^J(fre\\(frmmre\\,Ze),{rh\\,Zh)]x 
g-i(kiTe||-(-k2Tfc||)^ (A-1) 

where A denotes the QW surface area. Transforming 
to center-of-mass (CM) and relative coordinates in the 



QW plane (Rn — ar^u +/3r,j||,r|| — r, 



fell 



17l||: 



where 



a — nic/M^, (3 — 1 — a, and mc, Afx are the electron 
and exciton in-plane effective masses, respectively) we 
can decompose the exciton wave function into an enve- 
lope function, 4){r\\ , z^, Zh), and a free motion part related 
to the in-plane CM coordinate. Denoting kx as the in- 
plane CM momentum we find 



kl,k: 



XZe,Zh) 



Zh)e 



i[R|l ■(ki+k2)+r|| -(/Ski-Qka 



ikx-R 



II X 



e 

<5(kx 



ki - k2 



3ki-ak: 



Xz,,Zh) (A-2) 



We use the simplest exciton wave function for the exci- 
ton ground state: (l){r\\, Ze, Zh) = Mxe{ze)Xh[zh)e~'^''\^^ , 
where A/" is a normalization factor and A is a variational 
parameter associated with the exciton Bohr radius in the 
QW, which-is fixed by maximizing the binding energy of 
the excitorO. We note that the use of a wave function 
separable in z and r|| , facilitates considerably the calcu- 
lation of the scattering matrix elements, although it is 
strictly justifiable only for narrow well structures. As- 
suming perfect confinement of electrons and holes in the 
QW, and taking the z axis origin in the center of the 
QW, the confinement functions are: 

This approximation is quite accurate for almost all prac- 
tical thicknesses L. Calculating the binding energy of 
the IS exciton in a GaAs QW of width L — 200A, we 
find Eb = 6.73meV, A = 150A(£;b = 7.49meV, A = 133A) 
for the heavy-hole (light-hole) exciton. These results, ob- 
tained with our single parameter trial function, are quite 



close to those obtained inc3, which took into account fi- 
nite QW barriers, and anisotropic masses. 



Combining the first terms from Eqs. ([1.4), the direct 
term reads 



Vdir(g) 



32A2e2 /-^/^ r 

— —5- / dzedzhdzc / d kie ■= x 

[5(A|/3q - kil) - 5(A|aq + ki|)] ^(Afci) x 
cos^(7rZe/L) cos^i-KZh/L) cos^(7rZc/L), (A-4) 



where g{\q) was defined in Eq. (II. 6). We note that 
the electron wave functions in the QW plane are absent, 
as they contribute, together with the in-plane center of 
mass part of the exciton wave function, a fixed phase fac- 
tor which is unimportant for our purposes. The momen- 
tum integral is a simple convolution and can be readily 
evaluated 

j d^h [g(A|/5q - kil) - g(A|aq + ki|)] ^(Afci) 

^ J d^P<f'{p/>^) [exp(ip • /3q) - exp(ip • aq)] 
\g{XPq/2) - g{Xaq/1)] (A-5) 



47r2 
2A2 



where 



-p/X 



(A-6) 



is the exciton spatial wave function in the QW plane. 
In order to evaluate the z integrals, we change to the 
coordinates 

,_1 

Z — Zq Z^ , Z — ^ (Zg -|- Zc) 

which implies the following change in the integration lim- 
its: 

L/2 ^0 i-L+2z' 

dze / dzc — / dz' / dz 

-L/2 J-L/2 J-L/2 J-(L+2z') 

L/2 i-L-2z' 

dz' I dz, 

J-{L~2z') 

and performing the integrations results in the direct 
term, given by Eq. ([1.5) 



Appendix B 

In this Appendix we present an analytical calculation 
of the angular part of the exchange integral. This calcu- 
lation, while quite complicated, is worth the trouble for 
two reasons. The first is, obviously, a major reduction in 
the time consumed by the numerical computation, mak- 
ing it feasible. As a subsidiary benefit, the partially an- 
alytical calculation removes some convergence problems 
encountered in the previously used numerical algorithm. 

Using dimensionless momenta (e.g., by writing q we 
mean Xq) the exchange integral reads 
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/exc(Afc,g,6') = -An J dkikih{ki)d(f> \^[l + {ki- Ak + qf] [l + (q- Ak)2] ^^^j 



[l + (ki-Ak + aq)2] 



'3/2 



-At: / dkikih{ki){{I) + {II)} . 

r 



(B-1) 



where 61 is the angle between q and Ak and (j) is the angle 
between ki and Ak. We shall work out separately the 
two parts of Eq. (B-1). 
The first part is given by 



27r 



(/) = I d(f> [ad + [bd + ae) cos 4>— 
Jo 

c{d + a/ a) sini/) — c(e + bja) sine/) cos 4>-\- 

3/2 



he cos + c /a sin' 
where we have defined 



<f\' 



(B-2) 



a = l + kl + Afc^ + a^g^ - 2aqAk cos 6* 
5 = 2ki{aqcos9 - Ak) 

c — 2aqkismd (B-3) 



1 + A:f + Afc^ + r - 2qAk cos 6* 



e — 2ki{qcosd — Ak). 



Transforming 
we have 



(I) 



TT and then to z = tan((/)/2) 



dz{l + zy[{a-b){d-e)+ 

2c(d - e + a/a - b/a)z + 2{ad - be + 2c^/a)z^ + 
2c{d + e + a/a + b/a)z^ + (a + b){d + e)z**] "^^^ . 

Finding the roots of this 4th order polynomial we have 

/oo 
dz(l + z2)2 [(z2 + 2Az + B)x 
-oo 



{z"^ + 2Cz + D){a + b){d + e)] 



-3/2 



(B-4) 



where 



A 



5 = ^ ; C 



D 



a + b a + b a{d + e) d + e 

We note for later use that A, B, C, D obey the inequalities 

B>0 ; D>0 ; B - A^ > ; i:>-C^>0. 

(B-5) 



One can reduce the polynomial in Eq. ( B-4 ) to a product 
of two quadratic binomial forms, by making use of the 
following transformation 

p + qy 

z = 

1 + 2/ 

with p, q given by 

ID- B ll (D- B 
^''^=2:43c^VilA3c 



AD-BC 
A-C 



(note that we assume A^ C, deferring for later the treat- 
ment in this special case). Noticing that p and q always 
have opposite signs, we choose g < and divide the in- 
tegral into 



dz — 



dz + I dz 
'o 



where the limit z — oo corresponds to y = —\ and the 
limit z — —oo corresponds to y = —1 with q and p inter- 
changed. The integral now reads 



(/)^AA. 
where 



-1 [y^l + q^) + 2y{l+pq) + l+p^] 

dy— 

-p/q 



[(y2+s2)(y2+i2)]3/2 



(St)3 



2 _p^ + 2Ap + B 2 _P + 2Cp - 
* ~ g2 + 2Aq + B ' ~ q^ + 2Cq + D' 



1 r«/P^ [yHl+p^) + 2y{l+pq) + l + q^ 

1 ^ [(y2 + l/,2)(y2 + l/i2)]3/2 



D 



and 
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Making the transformation y = stanT/i in the first integral and y = \ tani/; in the second, we find 

1 2 



{I) 



St' 



tan ( 1/s) |^(gssin?/) +pcosV')^ + (ssin-)/) + COS-!/))^ 
"^"^ (1 - fc2 sin^ V')3/2 



tan ^ ( —p/sq) 

ta„-i(-?./p) r(p/ssinV' + gcosV')2 + (l/ssinV- + cosV')^] 
dip = 5 

tan-i(-s) (l-fc^sin ?/))3/2 



where we have defined 



fc2 = 1 - S^t^ 



(B-6) 



Changing -0 — > 7r/2 — ^ in the second integral, and using tan ^ s + cot ^ s = "n 12 we find that the second integral 
exactly matches the first except that its upper limit is larger by tt; thus we are left with 



{I) 



rr-tan s \ {I p'i) ^ir? ^ ^ {I + q^) cos^ ip + s{l + pq) sin2ilA' 
dip- — 

- tan~-^ s 



(1 - fc2 sin^ ^)3/2 

This integral can be evaluatedjising incomplete elliptic integrals (see p. 201 i Furthermore, one can exploit the 
functional relations (p. 911 iitB) 



F{^dp,k) = -F{4>) 
E{-cP,k) = -E{<p) 



F{mT -(p,k) = 2K(fc) - F{(P, k) 
E{n-K -<p,k)^ 2E(fc) - E{(P, k) 



where F{(p,k),E{(p,k) denote incomplete elliptic integrals of the first and second kind, respectively, and K(fc),E(fc) 
denote complete elliptic integrals of the first and second kind, respectively. The final result is 



(1) 



2N 



2{3 + 3p^q^ +4pq+p^ + q^ 



{s^ +r)F,{k) - 2t^K{k) 
1 



1 

2 I +2n 



(l+/)2 + (l + ,2)2^^2 

2s2E(fc) - +t2)K(fc) I . 

I 



(B-7) 



For i maginary values of fc, the result is the same as Eq. 
(B-7) with the i nterchanges s < — > t and k — > fc, k being 
defined in (B-6). 

We now return to e xam ine two special cases which can- 
not be treated using (B-7). The first case corresponds to 
A = C , B ^ D. Here we have 

p = zLco ; q = —A. 

Although s,t — oo, k is still finite, taking the value 



D^A^' 



thus, enabling the use of (B-7) with some modifications. 
Remembering the inequalities (B-5) we note that the sign 
of fc2 ig determined by the sign oi D — B. We obtain, for 
D> B 



A^D - A^ 



[(a + 6)(d + e)]^/' {D-B) 
(1 + ^2)2 



(B- 



[B-A^){D-A^) 



(\ ^ D - A^)¥.{k)- 



2K{k) - 2(3^2 + 1) 



2E(fc) - 1 



D-^2 



K(fc) 



For D < B the result is the same as Eq. ( B-8 ) with the 
interchanges B < — > D and k — > k. 

The second spec ial c ase corresponds to A — C ,B = D. 
Here the integral (B-4) reduces to 



dz- 



(1 + ^2)2 



{!) = 



This integral can be easily evaluated (see p. 81 inEil) and 
the result is 



[(a + 6)(d + e)]3/2 (z2 + 2Az + B)3- 



{!) = 



TT 3 + 4^2 + 25 + 3^2 

4 [{a + b){d + e)]'V2(B _ ^2)5/2 



(B-9) 



Finally, we evaluate the second part of the exchange 
integral, appearing in Eq. (B-1). We have 



{II) = 



1 



2ir 



(d— fc2)3/2yp (a + 6cos0 — csin0)3/2 ' 



where a, 5, c, d where defined in ( B-4 ) . Using the trans- 
formation (p — 2ip + where tan/3 — —c/b, we have 



(11) 



-13/2 



dip 



{d^klf/^J_p,2 (a+p-2psin20)3/2 
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where p = y/h^ + c^. This integral can be evaluated using 
the elliptic integral of second kind, resulting in 



{II) 



E 



2p 



Appendix C 



(B-10) 



Eq. ([II.l). It is useful to work with the Fourier trans- 
formed trion wave function (|III.6| ), where we have for the 
two spin configurations, explicitly 



In this Appendix we present the calculation of the trion 
binding energy, based on the trial wave function given in 



The effective mass Hamiltonian is given by 



2m, 



eo\/(ri - i"2)^ + (21 - ^2)^ 
(C-1) 



where a = rrie/ml, — nie/mf^ are the electron-hole effective mass ratios in the plane and z directions, respectively. 
Transforming to CM and relative coordinates in the QW plane: 



ftp 



cr(ri + r2) + r,i 



l + 2a 
rih = r,; - Th, i = 1,2 

and neglecting CM motion which does not contribute to the binding energy, we have 

2 

H ^fo(rih, Zi,Zh) + i?i2(rih, zi, r2h, ^2), 



where 



^fi2(i"ih,2;i,r2h,Z2) = -2o-V,.i,, • Vr^h + 



2 



\J (rih - r2h)^ + {zi - Z2Y 



(C-2a) 
(C-2b) 

(C-3) 

(C-4a) 
(C-4b) 



and we employed natural units of length and energy, namely, bulk effective Bohr radius and Rydberg. The expectation 
value of the trion energy is given by 



E : 



WW 



1±k2 



(C-5) 



where k was defined in Eq. ( III.4| ), and the upper (lower) signs apply to the singlet (triplet) spin configuration. The 
various terms appearing in Eq. ( |C-5| ) are calculated in the in-plane Fourier space, resulting in simple 2D convolutions, 
similar to those in appendix A. The results are: 



I^ = 



j^j [{<l>XiXh\Ho\<l>XiXh) + {(l>'xiXh\Ho\(l)'xiXh)] = ("2 + a,) + 



-{l + a)-X^ ( - 



dqqh{\q)giXq/2) 



+ [A ^ A'] 



I2 

Ji 
J2 



- {(^XiXh\HoWxiXh) ^ ^(t) (2 + ^.) + 



2 \L 



4(1 + g) 
(A + A')2 



^1X102X2|^12|01X102X2> = 2A3 



L 



dqqh{2\q)g{\q/2)g{\'q/2) 



6iXi0k2|i^l2|<^'iXl</'2X2) = 2/t2A3 / dqqh[2Xq)g\Xq). 



(C-6a) 

dqqh{\q)g{\q){C-&h) 
(C-6c) 
(C-6d) 
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In Eqs. (C-6), h{\q), g{\q) are the functions defined in Eq. (II. 6), and we have denoted A = A AV (A + A'). We note 
that in the 2D hmit (L — > 0), eliminating the kinetic energy term in the z direction, Eqs. ( |C-(]| ) reduce to 



Ii = 

Ji - 
J2 = 



1 

A2 



1 

A^ 



(l + .)-4(l + l 



4(1 + cr) 




8 


(A + A')2 


A 


+ A' 


( Ak 






(A'-A)- 




< A'fc 


11: 


+ h 


1 (A'-A)- 
^. 4A' 




Sttk 






2(A + A') 







E(l - AVA'2) - 2K(1 - AVA'2) 
E(l - A'VA^) - 2K(1 - A'VA2) 



A < A' 

A > A' 
A = A' 



(C-7a) 
(C-7b) 

(C-7c) 
(C-7d) 



where K, E are the complete elliptic integrals of the first 
and second kind, respectively. These closed form expres- 
sions coincide, with the integral expressions of Sandler 
and Proettocj, for a ^ and vanishing magnetic field 
(note a misprint in Eq. A6 of this paper). 

In order to maximize the trion binding energy, we must 
first write the expectation value of the exciton energy: 



1 

~o 
4A| 

L2 



(1 



dqqh{\)q)g{\^q/2), 



(C-8) 



resulting in the IS exciton Bohr radius and binding en- 
ergy values stated in Appendix A. Note that in the 2D 
limit we have 



1 



^2 



4 



(C-9) 



Minimizing Eq. (C-9) in the cr — > limit, results in Aq = 
aB/2, — — 4Ry, as expected. The trion binding energy 
with respect to the exciton is related to the trion and 
exciton energies through 



El' = E- 



En — Er, 



(C-10) 



where = (Tr/L)^ is the confined electron energy, in our 
infinite barrier model. We note that the kinetic energy 
terms in the z direction cancel out in the expression for 
the trion binding energy, so they can be removed from 
the calculation. Maximizing Eq. ( C-1C1| ) with respect to 
the two variational parameters A, A', we finally find the 
results given in section [II , where we used the GaAs bulk 
values of the effective Bohr radius (as = lOOA), and 
Rydberg (5.7meV). 
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